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It is shown, that a non-local form of the Gross-Pitaevskii equation 
allows to describe not only the long-wave excitations, but also the 
short-wave ones in the systems with Bose-condensate. At given 
parameter values, the excitation spectrum mimics the Landau 
spectrum of quasi-particle excitations in superfiuid Helium with 
roton minimum. The excitation wavelength, at which the roton 
minimum exists, is close to the inter-particle interaction range. 
It is shown, that the existence domain of the spectrum with a 
roton minimum is reduced, if one accounts for an inter-particle 
attraction. 



spectrum mimics the Landau one for the quasi-particles 
in the superfiuid Helium, which has the roton minimum. 
We study the influence of the inter-particle interaction 
on the dispersion curve in the short-wave domain. It is 
shown that the inter-particle attraction narrows down 
the existence domain of a stable spectrum. In particular, 
it reduces the range of parameters where a spectrum wi- 
th roton minimum may exist. 



The Gross-Pitaevskii equation (GP) obtained in pj [2j [3] 
is widely used in researches of superfiuid Bose-systems of 
weakly interacting particles. This equation describes the 
particle dynamics in a Bose-condensate. In particular, 
it gives a framework to study the condensate oscillati- 
ons and to obtain the elementary excitation spectra. The 
knowledge of the latter is necessary in order to obtain 
the thermodynamic functions of the system. 

At a stage of derivation of the GP equation for the 
condensate wave function $ (r, t) , it is usually assumed 
that <£> varies slowly over the distances of order ro, whi- 
ch is the range of the inter-atomic interaction potential 
U (r). This condition allows to express the GP equation 
in form of a non-linear differential equation. 

Unfortunately, in such a case, the imposed assumpti- 
on restricts a consideration to the long-wave excitations 
in the Bose-condensate. Namely, it allows one to derive 
the Bogolyubov sound (long-wave) oscillation spectrum; 
nevertheless, a (physical) superfiuid system can also have 
some well-defined long-lived short-wave excitations with 
the de Broglie wavelengths of the order of inter-particle 
distance (or interaction potential range). E.g., in the 
superfiuid helium the rotons are the excitations of this 
type. 

The purpose of this paper is to demonstrate that the 
non-local form of the GP equation describes the short- 
wave excitations in a system with Bose-condensate, as 
well as the long- wave ones, if the constraint on the 
of macroscopic wave function spacial change scale is 
released. At given values of the parameters the excitation 



NON-LOCAL FORM OF GROSS-PITAEVSKII 
EQUATION 

A system of many Bose-particles in the second quanti- 
zation picture with an account for the pairwise interacti- 
ons is described by the Hamiltonian 

t-2 



dri> + (r) ( -— A) V(r) 
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The equation of motion for the field operator in the Hei- 
senberg picture reads 

fc2 
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dip M) 
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Due to a large number of particles in the condensate, 
one may consider the operator ^(r,t) in eq. as a 
non-operator function ip{r,t). Further, it is convenient 
to define a function $ (r, t) : 



ijj (r, t) = $ (r, t) exp (— ifi t) . 
Then the equation reads 



d$ (ryt) 
dt 



2m" 



A + fx $ (r, t) 



(3) 
(4) 



+$(r,t) / dr'U(r -r')\$(r',t)\ , 
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where u is the chemical potential which value is determi- 
ned by the equilibrium of the system 



M = / dv'U (r-r')\$ (r')\' 



(5) 



with $o (r) the equilibrium value of the condensate 
wave function. Taking into account that the equilibri- 
um condensate density in absence of external fields no = 
|4>o(r)| 2 does not depend on the coordinate, one can 
express eq. (0} as 



d<$> (r,t) 



^— A 

2m 



n U $(r,i) 



(6) 



-$(r,i) J dr'U(r-r') |$(r',t)| 5 



where Uq — J dr' U (r — r'). If one is interested only 
in the long-wave perturbations, and assumes that the 
characteristic spatial scale of the condensate wave functi- 
on is much greater than the inter-particle interaction 
potential range, then the square of the wave function 
can be taken out of the integral in the right hand side of 
eq. (|6|). Thus, the GP equation takes the form 



ih 



<9$ (r,t) 

at 



— A + n Q U *(r,t) 
2m 



+U <Z>(r,t) |$(r,t)| 



(7) 



This is the usual way to apply the GP equation. Small 
oscillations described by eq. ([7]) have the dispersion law 



U n k 2 



h 2 ^ 

Am 2 



(8) 



where w, k are the frequency and the wave number. The 
dispersion law ((8]) was derived by Bogolyubov (4] in a 
different approach with the use of the second quantizati- 
on method. The excitation energy due to eq. ((8j) linearly 
depends on the wave number for small k and is monotone 
increasing with k, approaching the free particle dispersi- 
on law. 

One should stress that spectrum (j8]) is valid only for 
the long- wave oscillations due to the approximation used 
for derivation of eq. ([7]) . 

Meanwhile, the derivation of the GP equation in 
form (jnj) does not assume that the excitations are of 
long-wave type. Thus, eq. © is also valid for descri- 
ption of the excitations which wave length is of order of 
the inter-particle interaction. We call further eq. j6| the 
non-local Gross-Pitaevskii equation (NGP). This equati- 
on can be applied in order to study the dispersion law 
of the collective excitations of a Bose-condensate in a 



short-wave domain of spectrum and to analyze the spati- 
al behavior of the condensate wave function at the di- 
stances of order of the inter-atomic potential range. 



HYDRODYNAMIC FORM OF NON-LOCAL 
GROSS-PITAEVSKII EQUATION 



Equation (|6|) can be rewritten in the form of 
hydrodynamic equations, if the condensate wave functi- 
on is taken as $ = rje lip , where tp, rj are the phase and 
the absolute value of the condensate wave function and 
1 — n is the particle number density in the condensate. 



V 

On 
~dt 



Taking the above into account, eq. (j6J) leads to 



V(nv) 

h dip 
m dt 



0. 
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where the notion of superfluid velocity is introduced as 

v = (h/m)V<p. (11) 

Eq. ([9]) is the continuity equation, and eq. (flQ|) is 
analogous to the Josephson equation in the theory of 
superconductivity. Taking a gradient of both sides of 
eq. (JTDJ) and taking into account the definition for the 
velocity, one gets 



9v 

at + (vV)v 

where 



Vto, 



(12) 



-If dr 'u ( r _ r ') n ( r >\ 
2m rj m J 



In a linear approximation, the right hand side of 
eq. (|12j) can be expressed in form Vw = — Vp/mno, 
where the pressure and the density are related in a non- 
local way 



2 



?7oA?7 + n J dr'U (r — r') n(r'), 



(13) 



Non-local relation between the pressure and the 
density was used in the phenomenological approach of [S] 
(neglecting the quantum term) in order to describe the 
rotons in the superfluid Helium. From our considerati- 
on it follows that the non-locality kernel, which relates 
the pressure and the density in a phenomenological 
approach, is determined by the inter-particle interaction 
potential. 
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CO/COn 



SMALL OSCILLATION SPECTRUM WITH 
NON-LOCAL EFFECTS 

The non-local GP equation leads to the following di- 
spersion law for small oscillations 



U3 



Uknok h k 



2 7,4 



4m 2 



(14) 



which is different from ([8j due to appearance of 
a potential Fourier component Uk — J drU (r) e 1 kr 
(dependent on the wave vector) instead of Uq. 

Usually, at short distances, a strong repulsion occurs, 
thus the interaction potential grows rapidly. It leads to 
a divergence of the Fourier component Uk- In order to 
overcome this obstacle one may suppose that the potenti- 
al remains finite at short distances [6]. As the simplest 
case, we will use the „semitransparent sphere model" 
potential 



U(r) 



0, 



r < a, 
r > a, 



(15) 



and consider U m and a as two independent parameters 
of the dimension of energy and length, correspondingly. 
In this case the Bose-condensate excitation spectrum 
depends on the only dimensionless parameter 



rr 



8ir hq a 5 m U m ' 
and can be written in form 

^ 2 = w o/W, 



(16) 



(17) 




z=ka 



Fig.l. Condensate excitation spectrum for the „semitransparent 
sphere model" interaction; £ h = 3.4 ■ 10~ 3 , fj = 0.92 • 10~ 3 

At & < £ < £h the dispersion curve acquires a roton mi- 
nimum (middle curve in Fig.l) and its depth increases 
with decreasing £. At the value £ = = 0.92 ■ 10~ 3 , the 
excitation energy minimum equals to zero (lower curve in 
Fig.l); it means an instability of a many-particle system 
for given parameters of the interaction potential. Thus, 
in the range < £ < £/, the excitation spectrum mimics 
that of superfluid Helium with a roton minimum [7|. 

The roton minimum position z$ is determined by the 
equation /' (zq) — 0. Dependence of zq on the parameter 
£ is shown in Fig. 2. With increasing £ the roton mini- 
mum position gets shifted toward smaller wave numbers, 
changing from its maximum z max — 5.45 (for = 
0.92 -1CT 3 ) to minimum z mm = 3.72 (for £ h = 3.4- 10~ 3 ). 
It follows that the excitation wave length A, at the di- 
spersion curve minimum lies in the range 



where 

f(z) = z _1 (sinz — zcosz) + £z 4 , (18) 

ojq = An no a U m /m (19) 

z = ka. (20) 

Notice that the dimensionless parameter £ (j 1 6[) strongly 
depends on the potential range £ oc a~ 5 , rapidly growing 
with decreasing a. 

The dependence of condensate excitation frequency 
on the wave number is shown in Fig.l. At parameter 
values £ > ^ = 3.4 • 10 -3 , the dispersion curve does not 
have local minima and the energy monotone increases 
with wave number (upper curve in Fig.l). 



1, 15a < A < 1,69a. (21) 

Thus one may conclude that the existence of a dispersion 
curve minimum is caused by the finiteness of the inter- 
particle interaction range and is not connected with any 
specific superfluid properties of system. It is noteworthy 
that in Helium the repulsive range of the inter-atomic 
interaction potential may be estimated as a = 2.56 A 
and the roton wave length is A r = 3.3 A such that the 
relation A r = 1.28a fits into the interval (|2ip . 

Dependence of the roton gap A = Wny / (zo) of the 
spectrum (fTT)) . (fT8")) on the parameter £ is monotone; 
Fig. 3 shows that A /ujq is zero at £ = and achieves its 
maximum A max /uJo = 1.17 at £ = 
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A/CO, 




Fig. 2. Roton minimum position zq as a function of parameter £ pjg 3 Roton gap versus parameter £ 



ATTRACTION INFLUENCE ON THE 
SMALL OSCILLATION SPECTRUM 

In the majority of the papers devoted to a non-ideal 
Bose-gas, as in [4], it is supposed that the inter-particle 
interaction is of repulsive type. It is of interest to investi- 
gate the influence of the inter-particle long-distance 
attraction on the elementary excitation spectrum. 

Toward this end we calculate the spectrum for the 
interaction potential 



r (1-1 - > Umi r - a ' 

-C/r\ r>a, 



(22) 



where C is positive. This potential accounts for a Van 
der Waals long-distance attraction of particles. The di- 
spersion law for potential (f2"2l can be expressed as 



(23) 



with a dimensionless parameter a = C/U m a 6 and an 
auxiliary function 



z 2 

1 

6 



smz 



Z - 1- Z - 
3 V 2 



(24) 



Ml-** 



(Si (z) stands for integral sine) describes the impact of 
the attractive part of the inter-particle potential on the 
elementary excitation spectrum. This function is shown 
in Fig. 4. Function g(z) is positive for wave numbers 
< z < 2.45 and negative for 2.45 < z < 5.40. For 
z > 5.40 function g (z) is positive and rises sharply. 



One may notice that the interval 2.45 < z < 5.40, 
where function (|24|) is negative, overlaps with the domain 
of the roton minimum existence. It follows that the 
attraction leads to the growth of the excitation energy in 
the domain of the roton minimum existence and smears 
this minimum. In the range of small wave numbers, and 
also for very large z > 5.40, attraction leads to loweri- 
ng the excitation energy. Thus the long-distance atomic 
attraction results in smoothing of the dispersion curve 
and in narrowing the range of the existence of excitati- 
on spectrum, which has extrema (rotons and maxons). 
Attraction reduces the upper bound £/, and increases 
the lower one narrowing the range of the existence 
of a Landau like excitation spectrum <5£ = £/, — 
Dependence of 5£ for a spectrum with roton minimum 
on attraction strength a — C /U m a 6 , is shown in Fig. 5. 

Notice that experimental spectrum of quasi-particle 
excitations in the superfluid Helium at large momentum 
(p > 3A _1 ) increases slowly such that the derivative of 
co(p) is close to zero at the boundary of the spectrum [8]. 
This behavior can be understood with help of a long- 
range attraction. Recall that function g(z) is positive and 
increasing for z > 5,4. Thereby, due to (|23|) . it opposes 
the dispersion curve growing at large momentum. 

One should draw attention to smallness of the 
attraction part of the potential for Helium — it is the 
smallest among the rare gases. Namely, the depth of te 
potential well for Helium e = — 10.5K, for Neon e = 
-35.9K, for Argon e = -121K, for Krypton e = -173K. 
It seems that due to this fact the maxon-roton part 
of the dispersion curve in liquid Helium is the most 
pronounced, not only in a superfluid phase, but in a 
normal phase as well. Notice a feature of the superfluid 
phase: the lifetime of the maxon-roton excitations are 
much larger than that in normal liquids. 
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z=ka 

4 

Fig.4. Function g (z) in the domain of z of physical interest. 
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Fig.5. Range <5£ as a function of the interaction strength a. 

CONCLUSIONS 

It is shown that the Gross-Pitaevskii equation with 
account for the non-local effects, caused by the fini- 
teness of the inter-particle interaction range, can descri- 
be the short-wave excitation spectrum, for which the 
wave length is of inter-particle interaction range order. 

For a Bose-system the spectrum for arbitrary wave 
number is obtained; it has a shape of the Landau 
spectrum of a superfluid Helium. We estimated the range 
of parameters where such a spectrum exists. Computati- 
on shows the roton minimum in the dispersion curve is 
caused by the finiteness of the inter-particle interaction 
range; the wave length of this minimum is close to the 
repulsion range of the interaction. The influence of the 
long-distance inter-particle attraction on the dispersion 
curve is studied. It is shown that it results in smoothing 
of the excitation spectrum and in narrowing the range of 
the existence of a maxon-roton type curve. It is noticed 
that the attraction part of the potential for the Helium 



atoms is small in comparison with that for other atoms; 
it seems that this fact induces a pronounced maxon- 
roton part of the liquid Helium dispersion curve. 
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KOPOTKOXBHJIbOBI 3By/I>KEHH5I B 
HEJIOKAJIbHIH MO/IEJII rPOCCA-niTAeBCbKOrO 

A.n.IeaiuuH, KD.M.IIoAyeKmoe 
P e 3 K) m e 

IloKa3aHO, mo HejiOKajibHa cbopivia piBHHHpm rpocca-IliTaeBCBKoro 
fl03BOjisie onucara He tmbkh flOBroxBHjiBOBi, ajie it kopotkoxbh- 
jitoBi 36yssceHHsi b CHCTemax 3 6o3e-KOHfleHcaTOM. ilpn iigbhhx 
3HaneHHHx napaivieTpiB, cneKTp 36yfl}KeHB 6o3e-KOH^eHcaTy Hara- 
^ye cneKTp JlaH^ay KBa3inacTHHKOBHx 36y r a,>KeHb y Ha^njiHHHO- 
My rejiii" 3 potohhhm MmiivryMOM. ^ZIpBjKHHa xBHjii 36y r a,>KeHb, npn 
HKifi icHye potohhhh MimMyM, 6jiH3BKa ^o pamyca Mi^KHacTHHKO- 
Bo'i BsaeMOfli'i. IToKa3aHO, mo ypaxyBaHHH Mi>KHacTHHKOBoro npn- 
TaraHHH 3By»cye o6jiacTb icHyBaHHH cneKTpa 3 potohhhm MimMy- 

MOM. 



KOPOTKOBOJIHOBBIE B03MyiLIEHH5I B 
HEJIOKAJIbHOH MOflEJIH TPOCCA- IIHTAEBCKOTO 

A . II. MeaiuuH, K). M. HoAysKmoe 
P e 3 10 m e 

IloKa3aHO, hto HejiOKajiBHan cnopMa ypaBHeHHH rpocca- 
IiHTaeBCKoro no3BOjiHeT onncaTB He tojibko fljiHHHOBOjiHO- 
BBie, ho h KopoTKOBOjiHOBtie B036yjKfleHHH b CHCTeiuax c 6o3e- 
KOHfleHcaTOM. LTpH onpeflejieHHBix 3HaneHH5ix napaMeTpoB, 
cneKTp B036y»cfleHHH 6o3e-KOHfleHcaTa HanoMiiHaer cneKTp JlaH- 
flay KBa3HnacTHHHBix B036y»cfleHHH b CBepxTeKyneM rejiHH c po- 

TOHHBIM MHHHMy MOM . /IjIHHa BOJIHBI B036y>KfleHHH, npH KOTOpOH 

cymecTByeT potohhbih MHmiMyM, 6jiH3Ka k pa^nycy MejKiacTH- 
HHoro B3aHMOfleiicTBH5i. IloKa3aHO, hto yneT MejKHacraHHoro 
npHTSjKGHiiH cy>KaeT o6jiacTB cymecTBOBaHHH cneKTpa c potoh- 

HBIM MHHHMy MOM. 
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